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Abstract 

The dynamics of two coupled spins- 1/2 interacting with a spin-bath via the quantum Heisenberg 
XY couphng is studied. The pair of central spins served as a quantum open subsystem are initially 
prepared in two types of states: the product states and the Bell states. The bath, which consists 
of N (in the thermodynamic limit N oo) mutually coupled spins-1/2, is in a thermal state at 
the beginning. By the Holstein-Primakoff transformation, the model can be treated effectively as 
two spin qubits embedded in a single mode cavity. The time-evolution of the purity, z-component 
summation and the concurrence of the central spins can be determined by a Laguerre polynomial 
scheme. It is found that (i) at a low temperature, the uncoupled subsystem in a product state can 
be entangled due to the interaction with bath, which is tested by the Peres-Horodecki separability; 
however, at a high temperature, the bath produces a stronger destroy effect on the purity and 
entanglement of the subsystem; (ii) when the coupling strength between the two central spins is 
large, they are protected strongly against the bath; (iii) when the interaction between the subsystem 
and the bath is strong, the collapse of the two spin qubits from their initial entangled state is fast. 
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I. INTRODUCTION 

Compared with other physical systems {ll, sohd-state devices, in particular, the ultra- 
small quantum dots with spin degrees of freedom embedded in nanostructured materials 
are more easily scaled up to large registers and they can be manipulated by energy bias 



and tunneling potentials yjl. ^ 



or quantum computation 



I, a 



aturally, the spin systems are very promising candidates 
sl- Inevitably, the spin qubits are open quantum systems 



a, lZ| subjected to the interactions with their environments. In a short time, the states 
of the qubits will relax into a set of "pointer states" in the Hilbert space [18;], which can 
be quantified by using the purity [9]; the entanglement between the spin qubits will also 
vanish. Yet quantum entanglement is the most intriguing feature of quantum composite 
system and the vital factor for quantum computation and quantum communication 
0, l^- These two disadvantages, so-called decoherence and disentanglement, will not 
be overcome until the modelling of the surrounding environment or bath of the spin systems. 

For solid state spin nano-devices, the quantum noise, causing decoherence and disentan- 
glement of the qubits system, mainly arises from the contribution of nuclear spins, which 
is usually regarded as a spin environment. Recently, there are some works devoted to the 
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behavior of central spins under a strong non-Markovian influence of a spin-bath 
Lucamarini et al. had used a perturbation method [3] and a mean-field approximation 
to study the temporal evolution of entanglement pertaining to two qubits interacting with 
a thermal bath. They found entangled states with an exponential decay of the quantum 
correlation at finite temperature. Hutton and Bose |l^ investigated a star network of spins, 
in which all spins interact exclusively and continuously with a central spin via Heisenberg 
XX couplings with equal strength. Their work was advanced by Hamdouni et al. |l^, who 
derived the exact reduced dynamics of a central two-qubit subsystem in the same bath 
configuration. And they also studied the entanglement evolution of the central system. 
Yuan et al. jl7| 

coupled spins in quantum Heisenberg XY 
of all the above works are very exciting. Yet all of their methods are of complex analytical 



developed a novel operator technique to obtain the dynamics of the two 

1^ spin bath with high symmetry. The results 



derivations. And in Ref. 17|, their analytical results are dependent on some particular 
initial states. The study of quantum dynamics from different initial states, such as Bell 
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states and product states, is a very interesting issue in this system. Thus we introduce 
a "half analytical and half numerical" method here to solve this kind of open quantum 
problem and show light on their features of dynamics from different initial states. Moreover, 
the numerical part of our method is beyond the Markovian approximation [19] due to the 
strong non-Markovian behavior of such a center-spins-spin-bath model. 



In this paper, we study an open two-spin-qubit system in a spin star-like configuration 



which is similar to the cases studied in Ref. 



16 



171]: the interaction among the bath-spins, 



between the two qubits and between the subsystem and the bath are all of the Heisenberg XY 
type. The present model involves the Heisenberg XY interaction that has broad applications 
for various quantum information processing systems, such as quantum dots, Cavity-QED, 



etc [20, 



21 



22 



231 ]. First, we use Holstein-Primakoff transformation to reduce the model to 



an effective Hamiltonian in the field of cavity quantum electrodynamics 2J]. Second, we 
apply a numerical simulation to determine the dynamics of the whole system and obtain the 
reduced dynamics of the two coupled spin qubits by tracing over the bath modes. During 
our numerical calculation, there are no approximations assumed and the initial state of the 
subsystem (a pair of central spins) can be chosen arbitrarily. We will give some results 
about the purity, the z-component summation and the concurrence of the center open spin 
subsystem in the thermodynamical limit. Additionally, we will show that the bath can 
lead to entanglement of initially unentangled qubits. The rest of this paper is organized 
as follows. In Sec. [Tll the model Hamiltonian and the operator transformation procedure 
is introduced. In Sec. IIIIl we explain the numerical techniques about the evolution of the 
reduced matrix for the subsystem; Detailed results and discussions are in Sec. [IVl The 
conclusion of our study is given in Sec. |Vl 



II. MODEL AND TRANSFORMATION 



Consider a two-spin-qubit subsystem symmetrically interacting with bath spins via a 
Heisenberg XY interaction: both the subsystem and the bath are composed of spin-1/2 



atoms. Each spin in the bath interacts with the two center ones 

n — 

strength, which is similar to the cases in the literature [151. 



17 
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jy the same coupling 



251 ]. The interactions 
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between bath spins are also of the XY type. The Hamiltonian for the total system is 



H: 



H = Hs + HsB + Hb- 
Hs = /^o('S'oi + Sq2) + (5'cJi5'o2 
90 



SB 



N 



'5'oi'S'o2) 5 

TV 



1=1 



i=l 



N 



(1) 

(2) 
(3) 

(4) 



Here, Hs and Hb describe the subs yste m and bath respectively, and Hsb is the interaction 



part in the whole Hamiltonian 17|, |18|, |26||. /iq represents the coupling constant between a 
locally applied external magnetic field along the z direction and the spin qubit subsystem. 
Q is the coupling constant between two qubit spins. Sq^^ and Sq^^ (2=1,2) are the spin-flip 
operators of the subsystem qubits, respectively, which are: 









(0 




: 








: 













(5) 



and Sl[ are the corresponding operators for the zth atom spin in the bath. The index 
i in the summation runs from 1 to A^, where N is the number of the bath spins. Qq is the 
coupling constant between the subsystem and the bath, and g is the one between any two 
spins in the bath. 



Substituting the collective angular momentum operators J± = Sf into Eq. [3l we 
rewrite the last two parts of the Hamiltonian as: 



Hsb 
Hb 



9o 



2j 



(6) 
(7) 



where j = N/2. After the Holstein-Primakoff transformation (It transforms the spin bath 
of infinity spins into an effective boson bath) |27l |. 



(8) 
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with [6,6"*"] = 1, the Hamiltonian, Eqs. [6] and [71 can be written as 

HsB = go {Sqi + Sq2) 
h+h 
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h+[l-—]h + 



1 - —bb^ 
N 



1 - 



b+b 



In the thermodynamic hmit (i.e. 

HsB 
Hb 



^ oo) at finite temperatures, we then have 

9o[{S^i + S^,)b+{So, + S,,) 
2gb+b. 



(9) 
(10) 

(11) 
(12) 



Although the whole Hamiltonian composed by Eqs. [2], [TT] and [I2] is similar to that of a 
Jaynes-Cumming model there is an explicit difference between the two models. The 
present Hamiltonian describes two coupled qubits interacting with a single-mode thermal 
bosonic bath field, so the analysis of the model is a nontrivial problem in cavity quantum 



electrodynamics 20 



2l| . We note here that due to the transition invariance of the bath 
spins in our model, it is effectively represented by a single collective environment pseudo-spin 
J in Eq. [HI After the Holstein-Primakoff transformation and in the thermodynamic limit, 
this collective environment pseudo-spin could be considered a single-mode bosonic thermal 
field. The effect of this single-mode environment on the dynamics of the two coupled qubits 
is interesting and nontrivial. In Sec. IIV| we will show some results, for example, the revival 
behavior of the reduced density matrix and the entanglement evolution of the two central 
spins. This could be used in real quantum information applications. 



III. NUMERICAL CALCULATION PROCEDURES 

The initial density matrix of the total system is assumed to be separable, i.e., p(0) = 
® Pb- The density matrix of the spin bath satisfies the Boltzmann distribution, that 
is Pb = e~^^^'^ /Z, where Z = Ti (^e"^^^^) is the partition function, and the Boltzmann 
constant fc^ has been set to 1 for sake of simplicity. The density matrix p{t) of the whole 
system can be derived by 

p{t) = exp{-iHt)p{0)exp{iHt), (13) 
p(0) = Ps(0)®Pb(0), (14) 
ps{0) = \m){m\- (15) 
5 



In order to find the density matrix p{t), we follow the method suggested by Tessieri and 



Wilkie 



32l |. The thermal bath state Pb(0) can be expanded with the eigenstates of the 



environment Hamiltonian Hb'- 



M 



Pb(0) = XI \^rn)uJm{4>m\, (16) 
m=l 

g — -Bm/T 

^^m = , (17) 

M 

Z = J2 e"'''"^^- (18) 

m=l 

Here \(f)m), ^ = 1,2,3,- ■■ , M, are the eigenstates of Hb and Em are the corresponding 
eigen energies. According to the form of Eq. [121 I0m) = and Em = 2(7m, thus Um = 
exp(— 2(yfm/T) and Z = 1/(1 — e'"^^^^). With this expansion, the density matrix p{t) can be 
written as: 

M 

p{t) = ^UOm\^m{t)){^m{t)\. (19) 



m=l 

Where 

\^m{t)) = expi-tHt)\^mm = t/(t)|^„(0)). (20) 

The initial state is 

\^mm = \m)H. _ 

The evolution operator U (t) can be evaluated by different methods. In Ref . p/Zj , they use a 
unique analytical operator technique, which is dependent on the special initial state. Here, 



we apply an efficient numerical algorithm based on polynomial schemes [28|, |29|, |30| into 



this problem. The method used in this calculation is the Laguerre polynomial expansion 



method we proposed in Ref. 28| , which is pretty well suited to many quantum systems and 
can give accurate result with a comparatively smaller computation load. More precisely, the 
evolution operator U (t) is expanded in terms of the Laguerre polynomial of the Hamiltonian 
as: 

a+l oo / ■ , \ A; 



-w^(t^)" E(r^)-«-) 



k=0 

where a distinguishes different types of Laguerre polynomials [31]]; k is the order of the 
Laguerre polynomial. In practice, the expansion has to be cut at some value of fcmax for 
the compromise of the numerical stability in the recurrence of the Laguerre polynomial and 
the speed of calculation, /cmax is optimized to be 20 in this study and the time step t is 



restricted to some value in order to get accurate results of the evolution operator. At every 
time step, the accuracy of the results will be confirmed by the test of the numerical stability 
— whether the trace of the density matrix is 1 with error less than 10~^^. In a longer time 
scale, the evolution can be achieved with more steps. The action of the time evolution 
operator is calculated by utilizing recurrence relations of the Laguerre polynomial. The 
efficiency o 



in Ref. 



this scheme 



28| is about 10 times as that of the Runge-Kutta algorithm used 



32| . When the states |\E'm(^)) are obtained, the density matrix can be found by 



performing the summation in Eq. (fT9|) . 



In principle we should consider every energy level of the single-mode bath field: M oo. 
But the contribution of the high energy states \m),m > mc {mc is a cutoff to the spin bath 
eigenstates) is found to be negligible due to their weights Um as long as the temperature is 
finite. That is to say, Eq. (fT9!) could be changed to the following form: 

mc 

p{t) = Y,^m\'^n^{t)){^m{t)\. (21) 

m=l 

Given the density matrix of the whole system p{t), we can find the reduced density matrix 
by a partial trace operation, which traces out the freedom degrees of the environment: 

ps{t) = Tib (p(t)) . (22) 

For our model, ps = is a density matrix of the open subsystem consists of two 

coupled central spins, which can be expressed by a 4 x 4 matrix in the subsystem Hilbert 
space spanned by the orthonormal vectors |00), |01), |10) and |11). The most general form 
of an initial pure state of the two-qubit system is 

\iP{0)) = a\00) + (3\11) + 7|01) + 5|10), (23) 
with + + + |5|2 = 1. (24) 



IV. NUMERICAL SIMULATION RESULTS AND DISCUSSIONS 

After we obtain the reduced density matrix, we can calculate any physical quantities of 
the subsystem. In the following we focus our attention on three important physical quan- 
tities of the subsystem which refiect the quantum entropy increase caused by decoherence, 
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(a)Purity (b)(5o^+^02> 

FIG. 1: Time evolution of (a) purity and (b) z-component oscillation for an initial two-qubit state 
of IV'(O)) = |11) for different values temperature: T = d.2g (solid curve), T = 2g (dashed curve) 
and T = lOg (dot dashed curve). Other parameters are fiQ = 2g, go = g, Q = Og. 



the population inversion and the entanglement degree of the subsystem state respectively. 
These quantities are (i) the time dependence of the state purity, i.e. P = Tr(p|), which 
characterizes the conservation of the purity of the subsystem jsl, Q]. If P = 1, the subsystem 
is a pure state; whereas P = 1/2" (considering there are n spins-1/2 in the subsystem), it 
is in a completely mixed state ps = 1/2"J (J is the identical matrix in 2" dimension), (ii) 
the ^-component of the central spins, i.e. {S^-^ + 8^2), which demonstrates the population 



probability of the system, (iii) the time-evolution of concurrence 33|, |3J] for the two central 



spins of the open subsystem. The concurrence o: 



their intra entanglement, which is defined as 33l |: 



the two spin- 1/2 system is an indicator of 



C = max{Ai — A2 — A3 — A4, 0}, 



(25) 



where Aj are the square roots of the eigenvalues of the product matrix ps{o-^ ®<J^)p*s{o-^ ®cr^) 
in decreasing order. 



Theoretically, our method can deal with time evolution of the subsystem from any initial 
state. Here we first discuss the case of |V^(0)) = |11), which means both of the two center 
spins are in their excited states. The time evolution of the purity at different temperatures 



is given in Fig. 1(a) as well as the expecting value of the z-component in Fig. 1(b) We 
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FIG. 2: Time evolution of the concurrence for an initial two-qubit state of \tp{0)) = |11) at low 
temperature: T = 0.2g. Other parameters are fiQ = 2g, gQ = g, Q = Og. 

find that, (i) at a very low temperature, both the quantities present a collapse and revival 
phenomenon, which is identical with the two photons resonance of two two-level atoms in a 



cavity. In Fig. 1(a), the state of the first peak point {gti = 1.424) along the solid line is 



psigti) = 0.1472|11)(11| + 0.1472|00)(00| + 0.7056|fi) 

■ — I 

where \B) = l/v^(|01) + |10)). By the Peres-Horodecki separability [35] test, and if we 
make a positive partial transposition (PPT) operation to the 1st center spin, we will get a 
new matrix: 

pl'{gti) = 0.1472|11)(11| +0.1472|00)(00| +0.3528(|10)(10| + |01)(01| + |11)(00| + |00)(11|), 

whose spectrum is {0.500,0.3528,0.3628,-0.2055}. It shows explicitly that the bath can 
entangle the subsystem spins, although they are separable initially. And the state of the 
second peak point {gt2 = 2.468) is 

ps{gt2) = 0.9169|11)(11| + 0.0574|00)(00| + 0.02574|5) 

which could be approximated as its start state |11)(11|. In fact, one can find the time 
evolution of the subsystem concurrence in the condition of |^/'(0)) = |11) and T = 0.2g in 
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(a)Purity (b)(5o^+^02> 
FIG. 3: Time evolution of (a) purity and (b) z-component oscillation for an initial two-qubit state 
of \ip{0)) = 1 11) for different values 0,: 0, = Og (solid curve), 0, = 2g (dashed curve) and = 5g 
(dot dashed curve). Other parameters are fiQ = 2g, go = g, T = 5g. 
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(a)Purity (b)(5o^+^02) 
FIG. 4: Time evolution of (a) purity and (b) z-component oscillation for an initial two-qubit state 
of \ip{0)) = 1 11) for different values go: go = 0.2g (dot dashed curve), go = 9 (dashed curve) and 
5o = 5(7 (solid curve). Other parameters are /io = 2g, = Og, T = 5g. 

Fig. [21 It is interesting that the bath can induce entanglement between the two central 
spins periodically even without the direct coupling between the two spins. The peak state 
(t2 = 4:.972g) (the sohd line in Fig. l(a))is 

ps{t2) = 0.9992|11)(11| + 0.0008|00)(00| + 0.00002|S)(5|, 
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(b)C(O 



FIG. 5: Time evolution of (a) von Neumann entropy and (b) Concurrence for an initial two-qubit 
state of IV'(O)) = l/v^(|01) + |10)) for different values temperature: T = 2g (solid curve), T = 5g 
(dashed curve) and T = 20g (dot dashed curve). Other parameters are = 2g, go = g, = Og. 
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FIG. 6: Time evolution of (a) von Neumann entropy and (b) Concurrence for an initial two-qubit 
state of IV'(O)) = 1/^(|01) + |10)) for different values Q: n = Og (solid curve), n = 2g (dashed 
curve) and Q = 5g (dot dashed curve). Other parameters are fiQ = 2g, go = g, T = 5g. 

which could be approximated as its initial state |11)(11|. Then the subsystem revives at this 
time point and the period of the revival is about 5g. (ii) with increasing temperature, their 
oscillation magnitudes are quickly damped by the thermal bath: For the purity at T = lOg, 
P — s> 0.3 means a much large derivation from the initial state (to see the dot dashed curve 



in Fig. 1(a)); For the z-component summation of the two spins, it means the degeneration 
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FIG. 7: Time evolution of (a) von Neumann entropy and (b) Concurrence for an initial two-qubit 
state of \ip{0)) = l/\/2(|01) + |10)) for different values go: go = 0.2g (dot dashed curve), go = 9 
(dashed curve) and go = 5g (solid curve). Other parameters are ^xq = 2g, 0, = Og, T = 5g. 
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gt 

FIG. 8: Time evolution of (a) von Neumann entropy and (b) Concurrence for an initial two-qubit 
state of \ip{0)) = l/\/2(|00) + |11)) for different values go'- go = 0.2g (dot dashed curve), go = 9 
(dashed curve) and go = 55(solid curve). Other parameters are fio = 25, ^ = Og, T = 5g. 



of their magnetization (to see the dot dashed curve in Fig. 1(b)). Then we hold the bath 



at a moderate temperature T = 5g to research the role of the direct coupling between the 



two central spins. In Fig. 3(a), with increasing f2, the correlation between the two spins is 



strengthened so that the leakage of the information in the open subsystem is reduced and 
the collapse speed of the subsystem state slows down and distinct revivals are observed. In 
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Fig. 3(b), the whole magnetic moment oscillates around a mean value {Sqi + 5*02) = 0.6 
when Q is as large as 5g. One can see that if the direct coupling between the two central 
spins is strong enough, the two qubits can keep their initial state from the influence of the 
bath. The effect of the coupling strength between the qubit subsystem spins and bath spins 
Qq can be found in Fig. |H At a large value Qq = 5g, the strong interaction with the bath 
will quickly push the pure state of the subsystem into a mixed state (to see the solid curve 



in Fig. 4(a)); on the contrary, at a small value go = 0.2g, the initial state can be conserved 



to a great extent, P > 0.68, as the dot dashed curve shows in Fig. 4(b) 



Then we show the results of the temporal evolution for the quantum entropy 
and concurrence from a most entangled state (one of the well-known Bell states) 
\ip{0)) = l/v^(|01) + |10)) under the influence of the bath. As is known, the entropy 
of either particle inside a Bell states is 1, which is certainly in a most mixed state. 



On one hand, in Fig. 5(a), at any finite temperature, the entropy will return 1 after 



some fluctuation; on the other hand, in Fig. 5(b) , the entanglement degree between the 



two coupled spins will quickly collapse to zero and make some fluctuations and revivals 
accidently. The magnitudes of these fluctuations decrease with increasing temperature in 
both figures. The damping speed of concurrence evidently increases with temperature in 



Fig. 5(b) and when T = 20g, it shows a sudden death in a much short time and a tiny 



revival in a long time scale. In Fig. 6(a), the effect of subsystem inner-coupling f2 on the 



entropy is so weak that the biggest fluctuation of that is less than 0.08 and when Q = 5g, 



this effect can be neglected. Yet, in Fig. 6(b) , the concurrence can be restored to some 
extent with a comparatively larger Q: when Q = 5g, there is no sudden death of the 
entanglement between the subsystem spins and the concurrence oscillates about C = 0.4. 
From the above results, we find that the von-Neumann entropy is not a very good measure 
of this subsystem (consisted of two or more atoms) state as concurrence. In Fig. U\ with 
decreasing interaction between the open subsystem and the bath go, the influence of the 



bath spins on the subsystem is reduced. As is depicted in Fig. 7(b) ai go = 0.2g, the 
dot dashed line oscillates around the C = 0.7. This is evident that if we can keep the 
coupling between the subsystem and the bath as weak as possible, then the entanglement 
of subsystem can be conserved as much as possible. 
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With the combination of the technique of Holstein-Primakoff operator and our polynomial 
numerical scheme, we have thoroughly discussed the quantum dynamics of two central spins 
in a spin star lattice within XY interaction from different kinds of initial states. In principle, 
all the physical quantities can be calculated as the above ones. For the product states, it is 
noticed that the bath can help to entangle the two separable spins; For the Bell states, it is 
valuable that go, the interaction between the subsystem and environment, can be adjusted 
to slower the collapse of the most entangled state. However, we also find that the behavior 
of the reduced density matrix for the open subsystem relies on different initial state although 
they are all of the Bell states. By the comparison of Fig. [7] with Fig. [HI it is shown: in the 
former figure, when Qq is small, the concurrence can be maintained for a long time; while in 
the latter one, when qq is enhanced, it is quickly damped. 



V. CONCLUSION 



We have studied the quantum dynamics of the purity, z-component summation and con- 
currence of two coupled spin qubits in a apin bath via Heisenberg XY interaction. A novel 
numerical polynomial scheme is used in the calculation for the reduced density matrix of the 
central qubits after the model Hamiltonian has been performed by the well-known Holstein- 
Primakoff transformation. The procedure avoids the difficulty in Ref. that the initial 
state is limited by the analytic derivation. The time evolution of different types of initial 
states, either product states or Bell ones, is obtained. Although the subsystem is initially 
prepared in the product state (|11)) , it turns out that the bath can induce entanglement 
on the subsystem spins by the Peres- Horodecki separability test on the state p'^^(gti). On 
the other hand, the effect of different types of coupling on the entanglement of the system is 
studied. Generally, the coupling to the environment reduces the initial state entanglement. 
When the interaction between the two spins is large, they are protected strongly against the 
environment. Thus, it is found that the dynamics of the subsystem depends sensitively on 
the initial state, the bath temperature, the inner-coupling between the two central spins and 
the interaction between the subsystem and the environment. Besides, our numerical scheme 
is simple and independent of the initial state, which can be easily applied to the studies of 
some other kinds of open quantum systems. 
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